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I present a perturbative calculation of the spectrum of the Faddeev-Popov operator in Coulomb gauge in three 
dimensions, and Landau gauge in two and three dimensions, with an ansatz for the gluon propagator as the 
non-perturbative input. The results show how the low-lying Faddeev-Popov eigenvalue spectrum is modified as 
the first Gribov horizon is approached, and how the spectra can differ in Coulomb and Landau gauges. 



I. INTRODUCTION 

The low-lying spectrum of the Faddeev-Popov (F-P) oper- 
ator, in Coulomb and covariant gauges, is a probe of the in- 
frared properties of non-abelian gauge theories. Confinement 
in Coulomb gauge, in particular, is rather directly related to 
the F-P spectrum. The color Coulomb potential, for example, 
involves a product of inverse F-P operators, and the Coulom- 
bic self-energy of an isolated color charge, which is infrared 
divergent in a confining theory, depends crucially on the den- 
sity of low-lying eigenvalues of the F-P operator, as discussed 
below. The connection to confinement is less apparent in co- 
variant gauges, although the density of near-zero F-P eigen- 
modes is expected to be relevant to the infrared behavior of 
the ghost propagator. 

Coulomb and Laudau gauges are defined on the lattice as 
the set of elements, on each gauge orbit, for which the quantity 

^[f^] = -LETr[t/^W] (1-1) 

X ! 

is stationary with respect to infinitesimal gauge transforma- 
tions. Here d denotes the number of space dimensions, in 
Coulomb gauge, and the number of spacetime dimensions in 
Landau gauge, which is a convention I will adopt from now 
on. In general, along any gauge orbit, there are many station- 
ary points, known as Gribov copies, and at these points the 
F-P determinant may be positive or negative. This indefinite 
sign is closely to related to Neuberger's Theorem iQl], which 
demonstrates that BRST quantization of any non-abelian lat- 
tice gauge theory is ill-defined at the non-perturbative level. 
The picture is that in summing over all copies on a gauge or- 
bit, the copies with a positive F-P determinant are exactly can- 
celled by the copies with a negative F-P determinant, and the 
functional integral vanishes. It is for this reason that a con- 
straint of some kind, imposed on the domain of functional in- 
tegration, is necessary. Ideally the range of functional integra- 
tion should be a subspace (such as the Fundamental Modular 
region) containing only a single gauge copy with positive F-P 
determinant per gauge orbit, but at a minimum the integration 
range should lie within the Gribov region. This is the region 
which consists of all gauge copies in which the non-trivial 
eigenvalues of the F-P operator are all positive; i.e. the Gribov 
copies which are local minima of R\U]. These are, in fact, the 
configurations obtained by standard lattice gauge-fixing algo- 
rithms. The Gribov region is completely bounded, and the 
Fundamental Modular region is partially bounded, by the first 



Gribov horizon, where the lowest non-trivial J^'-P eigenvalue 
vanishes. It has been argued by Zwanziger |2(] that the vol- 
ume of the Gribov region is concentrated close to the horizon, 
much as the volume of a sphere in a high dimensional Eu- 
clidean space is concentrated near the surface. Since the di- 
mensionality of the space of all lattice configurations is very 
high indeed, the values of observables obtained at the Gri- 
bov horizon should dominate the expectation value. It would 
then be interesting to understand exactly how proximity to the 
Gribov horizon affects the behavior of various observables, 
starting with the spectrum of the F-P operator. 

As a step in that direction, this article presents a pertur- 
bative calculation of the F-P spectra. Perturbation theory is 
not necessarily trustworthy when dealing with the low-lying 
eigenmodes, but something may still be learned from it. In 
particular, it would be interesting to see whether proximity 
to the Gribov horizon changes the behavior of the low-lying 
spectra already at the perturbative level, and whether that be- 
havior is different, for some reason, in Coulomb and Landau 
gauges. The calculation is carried out for Landau gauge in 
two and three spacetime dimensions, and Coulomb gauge in 
three spacetime dimensions, to avoid the complications asso- 
ciated with renormalization in four dimensions.' The prox- 
imity to the Gribov horizon is controlled by a mass parameter 
in the transverse gluon propagator, which is where the non- 
perturbative information enters. I use an ansatz for the gluon 
propagator, motivated by Gribov's expression fl, which al- 
lows for any desired power behavior in the infrared. 



II. F-P EIGENVALUES AND THE COULOMB 
SELF-ENERGY 

The Coulomb potential between a static quark-antiquark 
pair located at points x and y is given by the expression 

yc{\x-y\) = ({M-')f^{-Vl){M-')>^) (2.1) 



Yang-Mills theory in Coulomb gauge is trivial in two spacetime dimen- 
sions if the spacetime manifold is fiat and non-compact, and for that reason 
Coulomb gauge in D = 2 will not be considered here. Non-trivial features 
associated with the Coulomb gauge F-P operator do appear even in two 
dimensions, if the space direction is compactified to S', and this case has 
been thoroughly discussed by Reinhai'dt and Schleifenbaum in ref. Q]. 



2 



where Cr is the quadratic Casimir of quarks in color represen- 
tation r, is the dimension of the adjoint representation of 
the gauge group, and M is the Faddeev-Popov operator, which 
is 



M'^ = -8«"V^ + gf''A'^{x)di W{x- y) 



(2.2) 



in the continuum. If Vc(|jc — yl) is confining, then this can 
only be attributed to an infrared singular behavior of M^', 
which must be related somehow to the low-lying F-P eigen- 
value spectrum. 

The perturbative evaluation of the F-P spectrum starts with 
the free-field, = Q case, on a finite periodic lattice of exten- 
sion L. The eigenmodes of the corresponding F-P operator are 
simply the plane wave states 
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(2.3) 



with a logarithmically rising potential. The criterion that the 
infrared divergence in the self-energy is stronger than loga- 
rithmic is 



> for some e > (2.8) 



This condition involves the near-zero F-P eigenmodes, as well 
as the eigenvalues. However, assuming that the eigenvalue 
spectrum is non-degenerate (apart from some rather special 
cases involving symmetric gauge-field configurations), then at 
fixed > each X„ is associated with a unique (n,A), which 
in turn determines p. Then — p^(k) in the infinite volume 
limit and 



lim 



> Hm 



(2.9) 



The proof of this inequality is given in the Appendix. It fol- 
lows that a sufficient condition for Coulomb confinement is 



lim/Pi^ 



> for some e > 



(2.10) 



and the Xa are some set of orthonormal vectors spanning the 
£(4 -dimensional color space. The F-P eigenmodes and eigen- 
values {(j)nji(-*') An,A} at > are also indexed by (n,A), 
denoted for brevity by n = (n,A), and it is assumed that the 
eigenmodes and eigenvalues are continuous and differentiable 
functions of g^, which smoothly approach (12. 3t as g^ — > 0. To 
connect the eigenmode spectrum to the Coulomb self-energy, 
we begin with the expression 



and inserting the spectral representation 



n 



this becomes Js]] 



X2 



where p(A-) is the normalized eigenvalue density 



1 



p(X)= lim YbiX-Xn) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



'Nc is the number of colors, and W ~ L . In 2h-1 dimensions, 
the integral in ( 12.61 ) is logarithmically divergent at the A- ^ 
end of the integration even in an abelian theory, and this is 
because the Coulomb potential in an abelian theory confines 



m. THE APPROACH TO THE HORIZON 

It was stated above that numerical simulations find local 
minima of R, which means, strictly speaking, that all of the 
eigenvalues of the F-P operator are positive. This statement 
has to be qualified a little. Even apart from Gribov copies, the 
Coulomb and Landau gauge conditions do not entirely fix the 
gauge, because if C/u(x) satisfies the gauge condition, so does 
GJ/u(x)G^, where G G SU(N) is any position-independent 
group element. This is a remnant global gauge symmetry, and 
it implies that at any stationary point of R there must be flat 
directions along the gauge orbit corresponding to zero modes 
of the F-P operator. These are the trivial eigenmodes 



(3.1) 



The statement that the F-P determinant is positive in the Gri- 
bov region really refers to the determinant of the operator in 
the subspace orthogonal to these trivial zero modes. 

Outside the Gribov region, some of the non-trivial F-P 
eigenvalues become negative, which means that for config- 
urations which lie exactly on the Gribov horizon there must 
be at least one non-trivial zero eigenvalue. However, in an 
infinite volume, the converse is not necessarily true: we can- 
not deduce, just from the fact that the spectrum of non-trivial 
eigenvalues begins at zero, that the gauge field lies on the Gri- 
bov horizon. Even in an abelian theory, which has no Gribov 
horizon, the spectrum of the F-P operator — in an infinite 
volume begins at A, = 0. 

Let us begin with g^ = 0, i.e. a free-field theory, with the 
eigenvalues and eigenstates shown in ( 12.3b . In this free case 
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(a) Type I scenario. (b) Type II scenario. 



FIG. 1 : Two scenarios for the behavior of the FP eigenvalue spectrum for gauge field configurations near the first Gribov horizon. Just outside 
the Gribov region (rf// < 0), there is a small interval of negative eigenvalues, which shrinks to a single eigenvalue at the horizon (rf// = 0). In 
the Type I scenario, the interval of negative eigenvalues begins at p = 0; at the horizon the nontrivial zero-mode is at p = 0, and at small p the 
eigenvalues grow with a non-standard power Xp ~ p^+'. For configurations inside the Gribov region {dn > 0) the growth Xp ~ is quadratic. 
In the Type II scenario, for configurations just outside the Gribov region, the interval of negative eigenvalues does not include p = 0, and at 
the Gribov horizon the non-trivial zero mode is at |p| > 0. 



we have 

p(X)-X(^-2V2 , (^,|-V2|(^,)=X (3.2) 

so that with an uhraviolet regulator, the Coulomb self-energy 
in J + 1 dimensions is finite for all space dimensions d >3, 
and marginally divergent (divergent as log(L) as extension 
L — > oo) at = 2. The latter divergence is expected, since 
the Coulomb potential increases logarithmically in 2 + 1 di- 
mensions, so the question in 2h-1 dimensions is whether the 
condition in eq. (12.8b is satisfied for some e > 

Outside the Gribov region, some of the non-trivial F-P 
eigenvalues become negative, and approaching the first Gri- 
bov horizon from the outside, the range of negative eigenval- 
ues should shrink away. Right on the horizon there must exist 
a non-trivial zero eigenvalue even for a finite spacetime vol- 
ume. So let us imagine increasing away from zero, and also 
placing a constraint in the functional integral by introducing a 
dimensionful parameter dn, and requiring that if dn > 0, the 
integration is over gauge fields inside the Gribov region, lying 
a distance djj from the first Gribov horizon, while if dn < 0, 
the integration is over gauge fields outside the Gribov region, 
at a distance \dH\ from the horizon. Then 

(Xp) = Xf^ + {AXp) 

^ p'-{l~F[g,p,dH]) (3.3) 

We use p as an index because, in the infinite volume limit, it is 
better to replace the integer index n by the continuous index p. 



Also the expectation value of A-p ^ can depend on neither the 
index A, since this would violate global color symmetry, nor 
on the direction of p, which would violate rotation invariance. 



If ^ = then F = 0, but we may speculate on the behavior 
of F[g,p,dH] atg^ >0 as dn varies. Suppose F has the form, 
near p — 0, 

F[g,P,dH] = a[g,dH]-b[g,dH]p' 

+higher powers of /7 (3.4) 

andb[g,dH] is positive for small Ate/// > all non-trivial 
eigenvalues are positive, so it must be that a[g,dH] < 1 for 
small p. Note that the eigenvalue spectrum in an infinite vol- 
ume still starts at A, = 0, even though the configurations are, by 
definition, off the Gribov horizon. At J// < some eigenval- 
ues are negative, and if those are the eigenvalues near p = 0, 
it means that a[g,dH] > 1. The negative eigenvalues must just 
disappear at dn = 0, and this is obtained if a[^,0] = 1 exactly. 
In this last case the subleading power of p in F[g,p,dH] takes 
over, and we have 

K ^ p'+' 

p{X) ^ X^d-2-s)/i2+s) (35) 

This a qualitative change in the low-lying F-P spectrum, com- 
pared to the behavior inside the Gribov region, and the suf- 
ficient condition ( 12.101 ) for Coulomb confinement is satisfied 
if 



^ To compute p{X), begin with the volume measure in momentum space, pro- 
portional to p'^^^dp, and change variables to ^ = X(p) to aiTive at p(X)dX. 



2s + 2>d (3.6) 
Inside the Gribov region, at p — > 0, the spectrum is simply a 



4 



rescaling of the zeroth-order spectrum 

Xp = {\-a[g,dH\)p^ 



(3.7) 



and, in the case of Coulomb gauge, the Coulomb self-energy 
is finite. The conjectured behavior of (kp) vs. p, for cIh posi- 
tive, negative, and zero, is sketched in Fig. |l(a)| 

But the scenario just outlined is not the only possible be- 
havior near the horizon. Consider, in particular, the case that 
b[g,dH\ is negative for small |c///|. Then we have 



- higher powers of p 
(3.8) 

and it is possible that (X^) is positive near p — Q) where the p^ 
term dominates, but negative in some finite region away from 
p = Q. The conjectured behavior in this case, for positive, 
negative, and vanishing du, is indicated in Fig. |l(b)[ and in 
this case we would still have A-p p^ at the horizon, for small 



where 

A'; - ^TT[aaiUj{x) ~ U]ix))] (4.2) 

The dimensionless lattice coupling gi is related to the gauge 
coupling ghy gj^ — a^^^g^, where a is the lattice spacing and 
D is the spacetime dimension. The eigenvalues and eigen- 
vectors of Kq are those shown in eq. ( 12.31 ). The operator M\ 
vanishes in the continuum limit, so I will just ignore it in what 
follows, and treat K\ as the only perturbation to Kq. Lattice 
Fourier transforms will be defined symmetrically 



-ikx 



(4.3) 



Of course, quantization in Coulomb and Landau gauge does 
not involve setting du to some definite value. What is re- 
quired, however, is a constraint on the range of functional in- 
tegration to lie within the first Gribov horizon. If it is true that 
entropy dominates due to the high dimensionality of the con- 
figuration space, and almost all of the volume of the Gribov 
region is concentrated at the horizon, then only lattice config- 
urations at or very near the horizon will contribute to vacuum 
expectation values in Coulomb and Landau gauge, just as if 
the constraint dn =0 were imposed. 



The first-order correction to X.p''' is 



AX 



(1) 



- {pA\Ki\pA) 



1 r 



1 



.■^X^e'"'^xi^LAf(0)sin(p,) 



(4.4) 



IV. PERTURBATIVE EVALUATION OF THE F-P 
SPECTRUM 



Now, according to the above definition of the lattice Fourier 
transform, the lattice A-field at zero momentum is 



The possible spectra shown in Figs. |l(a)| and |l(b)| are pure 
speculation at this point, but it is interesting, and somewhat in 
the spirit of Gribov's original work to see how far we can 
go in understanding the F-P spectrum with ordinary perturba- 
tion theory. 

Let us begin with lattice SU(2) gauge theory in either d+1 
spacetime dimensions (Coulomb gauge) or d spacetime di- 
mensions (Landau gauge), starting on a finite t/-dimensional 
volume V and taking the infinite volume V ^ °° and lattice 
spacing a ^ limits at the end. The F-P operator on the lat- 
tice is given by Jst] 



< 
(Kit 



{M,r^ = -8«'^E{8vv[(l-iTrf/,(x)) + (l-iTrf/,(x-/))] 
-5.+?,,(l - 5Trf/,-(x)) - 8,_.,,(1 - iTrf/,(y))} 



(4.1) 



(4.5) 



with —2/g < Af (x) < 2/g. Then suppose that the lattice A- 
field in Coulomb or Landau gauge has a finite correlation 
length /. This implies 



(4.6) 



where J? is the average value of Aj' in a hypercubic region 
of volume f'. Then, because of the factor of 1 /VV in (I4.4l i. 
the first-order correction to A-^ vanishes in the infinite volume 
limit. Of course, the first-order contribution vanishes even 
in a finite volume upon taking the expectation value, since 
(Af(0))=0. 



At second order 



AAp.A = 



L 

k.B 



\{k,B\Ki\p,A)\"^ 



(4.7) 
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where 



form of the full, equal-times gluon propagator. This gives 



ik,B\Ki\p,A) 



1 



(4.8) 



Then 



1 1 



:Y^A'l{k-p)A]{p-k) 



+ e-'^') (-e-'Pi + e''-') (4.9) 



In preparation for taking the continuum limit, we need to 
indicate powers of the lattice spacing explicitly. Let 



AX = a AX' , p^ap' , A' (x) = aA'f (x 



(4.10) 



where the primed quantities have the standard engineering di- 
mensions of these quantities in the continuum formulation. 
We also have, using Ak' = 2k/ (La), 



Ak'V 1 



mLD.jip'-k') 



2n J I'W _ I'io) 
X (^-e'P'i" + e-"''i"^ (-e-'P'f + e''^^') (4. 14) 

At this point we can take the continuum limit, and make use 
of the transversality property qiDij{q) = of the gluon prop- 
agator, to obtain^ 

(Ax;,) = g^Y^{xl^"''x%){xi^"'h{) 

B 

The primes, having served their purpose, will now be 
dropped. It is understood that the unprimed quantities now 
have their standard engineering dimensions. 

Using the competeness property 

IXbX^ = 5^'^ (4.16) 

B 

we sum over the color indices, which just gives an overall 
factor of two. The result is 



1 1 



Ak' 
2n 



Inserting these identities into (|4.9l l 



- y^lL{xl^"''X%){x'i^''h{) 



Ak'V 1 



271/ „2n'(0)_l'(0)^ 



(4.11) 



l^A>^{k~p)A^j{p~k) 

•J 

X (^-e'P'-" + e-'^'i"^ (^-e-'P'f+e"'''") (4.12) 



iAX„)^~2gpu,,J j-^-^^^„ 



D,jip-k) (4.17) 



(note the interchange of k^ and p^ in the denominator). 
Changing variables to q = p — k, and writing 



D,{q)={5,j-^ 



D{q) 



gives 



{AXp) = -2^2 



d'^q D{q) 



(27t)^ [p-qY-p^ 



(4.18) 



(4.19) 



We now go to (/-dimensional spherical coordinates 



Now we take the vacuum expectation value of AA,', noting 
that 

(Af (/t)A^(-fc)) = a^-''8'''Dij{k') (4. 13) 

where, in Landau gauge, Dij{k') is the full (i.e. dressed) gluon 
propagator In Coulomb gauge it is the spatial Fourier trans- 



d'^q = A,/_i / dq / sin''"^ q (4 20) 



We have ignored, in lattice regularization, tlie case that =4"", which 
would have to be handled by degenerate perturbation theory. This case is 
zero measure in the continuum limit, and will not require special treatment. 
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where 



and 



Define 



D{q)^q'-^D{q) 



(4.21) 



(4.22) 



Then 



2Ad-i 
{2ny 



l/(27t2) d: 

l/(67t3) d: 



-2 
-3 

:4 



(4.23) 



where 



(Alp) = -g^R,, Tdqq''-^ rdBsm''-^d{l-cos^d)^-^ ^D{q)p^ 

Jo Jo q^-2pqcosQ 

= -g^Rd r dQ &m'^-^-Q{\-cof?Q) [ dq J- -^q'^'^D{q)p'- 

Jo Jo ^-2pcos9 

= -g^RdpHh+h) (4.24) 



h = r'^dB sm''-^B(l-cos^B) f dq J- -D(q) 

Jo ^ 'Jo ^q-2pcosB 



h = I t/e sin'^^-efl -cos^e) / dq ^ D(q) 

A/2 Jo q-2pco^B 



Make the change of variables 9 ^ 7i — 9 in /2 

l2= J9 sin''"2 9(1 -008^9) / dq -D(q) 

Jo Jo ^ + 2pcos9 

Then, in Ii, it is useful to rewrite the integral over momenta q 

dl = / + / + / rdq -^^^ 



2pcos9^ D(^) ^ Z'^''™'*'^ D(4pcos9-^) _ D(4pcos9 + ^) 



(4.25) 



(4.26) 



^ — 2pcos9 

— 2pcos9 Jo 2pcos9 — (7 ^ 2pcos9 + ^ 

2/7COS9 — ^ Jo 2/:icos9 + ^ 

Altogether, we have to second order 

(Xp) = 4"' + (AXp) = p2(l -g'^R,I[p,mM) (4.28) 

where 

/■'t/2 f Z"" 1 ~ ~ 

/b,m,a] = / £/9 sin''"2 0n _£.os2 0N I / [D(4/7Cos9 + ^) 

Jo iJo ^ + 2/9cos9 

2pcos9 \ ^ ^1 

dq- [D{4pco&B-q)-D{q)]\ (4.29) 

2pcos9 — J 



I 

The m, a in I [p,m, a] are constants I will use to parametrize V. AN ANSATZ FOR THE GLUON PROPAGATOR 

the transverse gluon propagator D(^). 

The gluon propagators Dij in Coulomb and Landau gauges 
are transverse with respect to spatial momenta q in + 1 di- 
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mensions, and spacetime momenta in d Euclidean dimen- 
sions, respectively. Therefore these propagators have the form 
shown in ( 14.181 ). In a free theory 



of gluon fields 



D{ci) = 



1/(2(7) Coulomb gauge 
\/q^ Landau gauge 



(5.1) 



where the propagator in Coulomb gauge is at equal times, with 
q the space (rather than spacetime) momentum. The behavior 
( 15. lb is expected at high momenta, but it is certainly not cor- 
rect at low momenta, as seen from lattice Monte Carlo sim- 
ulations. In Landau gauge, the current evidence is that D{Q) 
is finite and non-zero at ^ = in three and four dimensions 
iBi, while ^0 in two dimensions M. In Coulomb 
gauge it appears that D{q) — > in four dimensions ifioll . 

In order to allow for non-singular power behavior in the 
transverse gluon propagator as p 0, 1 will adopt the ansatz 
that 



D{q) 



in Coulomb gauge, and 



1 



27? 



D{q) 



,2+a 



Iq^ 



(5.2) 



(5.3) 



in Landau gauge. Gribov's proposal for the gluon propagator 
in these cases corresponds to a = 2. The propagators go over 
to free-field behavior ■&& q^<x>. 
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FIG. 2: Equal-times Coulomb-gauge gluon propagator in 2+1 dimen- 
sions, at p = 6 and lattice volume, for L = 24, 32, 50. 



I am not aware of any lattice Monte Carlo computation of 
the transverse gluon propagator in Coulomb gauge in d = 3 
dimensions, in position space. In Fig.|2]l show data for D{R) 
obtained from the equal times correlator 



(Tr[A,(x,r)A,-(y,0] 



Aj{x,t)^-iUj{x,t)~u]ix,t)) 



(5.5) 



on the lattice. The correlator is calculated via lattice Monte 
Carlo with an SU(2) Wilson action, on an L? lattice volume at 
coupling p = 6 and L = 24, 32, 50, with the equal-times corre- 
lator computed after transforming the gauge fields to Coulomb 
gauge. Note that as the lattice volume increases, the gluon 
propagator develops a "dip" and actually becomes negative at 
the larger R values. This behavior appears to rule out a = 0, 
in which the propagator should be everywhere positive. A re- 
liable computation of D{q) as (7 — > will probably require a 
large-scale lattice calculation, as has been done for the Landau 
gauge. 



VI. RESULTS FOR F-P SPECTRA 

In section |lll] I introduced a parameter dn to control the 
approach to the first Gribov horizon, and speculated on the 
low-p behavior of Xp as the horizon is approached. In the 
perturbative calculation, the mass parameter m in the gluon 
propagator plays essentially the same role as dn- Note that 
in dimensions lower than 3h-1, where I[p,m,a] is convergent, 
the coupling is dimensionful, and we may as well choose 
units such that g~ = 1. Then 



{Xp)=p^{l-R,I[p,m,a] 



(6.1) 



Expanding I[p,m,a] in leading powers of p near p = 0, we 
have 

Rdl[p,in,a] = a[m,a] — b[m,a]p''^ 

+higher powers of p (6.2) 



in which case 



(1 - a[m,a])p^ + b[m,a\p^^^ 

+higher powers of p (6.3) 



(5.4) 



Suppose, for a given a, it is possible to find a critical value 
m = nic such that a[m£., a] = 1 ?indb[m, a] > 0. In that case we 
have the Type I scenario conjectured in Fig. |l(a)| above; i.e. 

1. m < nic and a[m,a] > 1; The low-lying F-P eigenvalue 
spectrum has a range of negative eigenvalues, starting 
at p = 0. We interpret this to mean that the trans- 
verse gluon propagator, which determines the spectrum 
at second order, is determined by configurations outside 
the Gribov region. 

2. m ~ nic and a[mc^a\ = 1: The region of negative eigen- 
values just disappears, and A,p ~ p'^^\ This is the case 
of particular interest, where the gluon propagator is de- 
rived from configurations which mainly lie right on the 
Gribov horizon. 
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Landau D=3 no solution Type I 

Landau D=2 ^o solution 

Coulomb D=3 no solution 



Type II 
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Type ^^^^^1 
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^ Type II J 
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FIG. 3: Summary of the qualitative behavior of the low-lying F-P spectra, according to 2nd order perturbation theory, for Landau and Coulomb 
gauges in D=2 and 3 dimensions. The sketch illustrates how the behavior of the F-P spectra depends on the assumed infrared behavior of the 
gluon propagator, which is parametrized by the exponent a. 



3. m > nic and a[mc,a] < 1. In this case the low-lying 
spectrum Xp — {1 — a[m,a])p^ is just a rescaling of the 
free-field spectrum, and the gluon propagator is derived 
from configurations inside the Gribov region. 



It should be noted at this point that the Type I scenario is 
in some ways reminiscent of the Dyson-Schwinger approach, 
and indeed eq. ( 14.28b resembles the Dyson-Schwinger equa- 
tion for the ghost propagator in covariant gauges (see, e.g., 
Fischer iflTll ). Of course these equations are not the same; 
(16.1b is an equation for the expectation value of F-P eigenval- 
ues, not the inverse ghost propagator, and it is derived from a 
perturbative expansion, not the Dyson-Schwinger equations. 
Nevertheless, the scaling solution |[l2ll is obtained from the 
Dyson-Schwinger equation by tuning a coupling so that the 
bare inverse ghost propagator in that equation is exactly can- 
celled by another term. In the absence of this tuning, the de- 
coupling solution ifisll is obtained. Similarly, in our approach, 
a mass parameter is tuned to exactly cancel the p-^ term in the 
eigenvalue spectrum, resulting in an enhanced density of near- 
zero eigenmodes. The motivation for the tuning in our case is 
to study the F-P spectrum at the Gribov horizon, which is only 
relevant to physics if, in fact, the functional integral over the 
Gribov region is dominated by horizon configurations. 

The Type II scenario is obtained if b[m, a] is negative when 
a[m,a] = 1, in which case there is still a range of negative 
eigenvalues, so this value of m is not the critical value. The 
critical value, corresponding to the horizon, is obtained at a 
value m = nic where fl[m,a] < 1, such that the function 



1 — fl[mt,(x] \p — 



b[mc,a] p' + c[mc,a]p'' (6.4) 



approximating (kp) at small p has a zero value, but no nega- 
tive values, for one choice of p ^ 0. In this case the horizon 
does not alter the power dependence Xp ^ p^ near p = 0. 



Both the Type I and Type II scenarios assume that 

a[m,a] = Rdl[p = 0,m,a] (6.5) 

is finite. This is not necessarily the case, however, and it is 
easy to check that l[0,m,a] is divergent at all a < for Lan- 
dau gauge in two dimensions and Coulomb gauge in three di- 
mensions, and is divergent for all a < ^ 1 for Landau gauge in 
three dimensions. There is no choice of m, for those choices 
of a, which completely eliminates negative F-P eigenvalues. 
This will be referred to as the "no solution" case. 

In order to determine which scenario is realized, at each 
choice of a for which /[0,m,a] is finite, it is necessary to cal- 
culate I[p,m,a] numerically. The result, for Coulomb gauge 
in three dimensions, and Landau gauge in two and three di- 
mensions, is indicated schematically in Fig. |3] To illustrate 
how these results are obtained, we consider in particular the 
case of a = 1 for Coulomb and Landau gauges in three di- 
mensions (i.e. d = 2 for Coulomb, and c/ = 3 for Landau). We 
begin with Coulomb gauge (Figs.|4]|7]i. Figure [4(a)] shows the 
low-lying F-P spectrum at a = 1 and m = 0.20 < nic, and it is 
clear that there is a region of negative eigenvalues starting at 
p ~Q. As m is increased, the region of negative eigenvalues 
shrinks in size, until at a critical value m = mt(oc) the interval 
of negative eigenvalues just vanishes. Figure [4(b)| displays the 
low-lying spectrum just below, at, and just above the critical 
mass at a = 1, which is nic = 0.2228. At m ^ nic, Xp is propor- 
tional to p^ near p = 0, with a proportionality constant which 
is positive or negative, depending on whether m is greater or 
less than nic- But precisely at m = nic, we find that Xp °^ p^^\ 
with s = s{a) > 0. Fig. |5]is a log-log plot of Xp vs. p over 
a large range of p, at a = 1 and nic = 0.223. For the range 
< /9 < 1, we can determine that s — 0.53 in this case, and 
Xp w 1.21/7^-^^ at small p. At around /? = |p| = 1 (in units 
g'^ = 1), the power behavior shifts to the free case, Xp = p^, 
and continues that way for all higher p, as expected. This is 
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FIG. 4: F-P spectra at a = 1. (a) m = 0.20 < nic- There is an interval of negative eigenvalues in the region < p < 0.009. (b) Xp at low p, for 
m = 0.20 < itic, in = nic = 0.2228, and m = 0.25 > nic- 




05 0.0001 0.001 0.01 0.1 1 10 100 



FIG. 5: Log-log plot of the spectrum of the Fadeev-Popov operator, 
for a = 1 at the critical = 0.223. A best fit at p < 1 yields A.^ = 
1.21p"3. 



an example of the Type I scenario. 

The next question is how ;7z^ and 2 + s change as a is varied. 
As already noted, we must choose a > to reach the hori- 
zon, which means that D{0) = 0, and therefore the transverse 
gluon propagator must vanish at zero momentum for Coulomb 
gauge in 2+1 dimensions, and for Landau gauge in 1h-1 di- 
mensions. As a — > 0+, the increasingly singular behavior of 
the integrand in /[p,m,a] must be countered by an increas- 
ingly large value of m^ , in order to satisy a [nic] , oc] = L A plot 
of nic vs. a is shown in Fig. |6] 

The power behavior Xp = bp^^" in the low-lying spectrum 
is crucial for Coulomb confinement, and the exponent 2 + s 
vs. a, obtained at ni = nic is shown in Fig. |7(a)| In 2+1 di- 
mensions the condition for Coulomb confinement (beyond the 
marginal divergence of the free theory) is that s > 0, which is 
seen to hold throughout the range shown. 
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FIG. 6: Critical value iric for the mass parameter in the transverse 
gluon propagator, vs. the power a. 



We also see that there is a sudden jump in s from roughly 
i = 1 to = 2 at a = 2. This is where the transition from Type 
I to Type II behavior takes place. As a 2 the coefficient 
b[mc,a] approaches zero (cf. Fig. |7(b)[ ) and then changes sign. 
Exactly at a = 2, where b[mc,a] = 0, the term which has the 
next higher power in p takes over, accounting for the sudden 
jump in s. 

Landau gauge in three dimensions, at a = 1, furnishes an 
example of the Type II scenario. The F-P spectrum at small p 
is shown in Fig.[8]for the mass parameter above (m ~ 0.087), 
below (m = 0.086), and equal m = mc~ 0.08644 to the critical 
value. 



VII. CONCLUSIONS 

If the integration over gauge fields is dominated by con- 
figurations on or near the first Gribov horizon, then the low- 
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FIG. 7: (a) Exponent 2 + s vs. a; and (b) coefficient b vs a; for the power-law behavior Xp - 
for the Coulomb gauge F-P spectrum in 2+1 dimensions. The sudden rise to 2 + s = 4 at a ^ 
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FIG. 8: The low-lying F-P eigenvalue spectra near the Gribov hori- 
zon, for Landau gauge in D=3 spacetime dimensions and a = 1, ac- 
cording to second-order perturbation theory. This is an example of 
the Type II scenario. 



est non-trivial F-P eigenvalue must be very close or equal to 
zero, even in a finite spacetime volume. The main finding of 
the perturbative treatment presented here is that if there is, in 
fact, a non-trivial zero mode, and the F-P eigenvalues are la- 
beled by the lattice momenta, then this non-trivial zero mode 
may occur at either zero momentum (Type I scenario) or non- 
zero momenta (Type II scenario), depending on the infrared 
behavior of the gluon propagator. While the spectrum of F- 
P eigenvalues does not translate directly into a prediction for 
the behavior of the ghost propagator (because the momentum 
behavior of the F-P eigenmodes must also be taken into ac- 
count), it is natural to conjecture that the Type I scenario is 
associated with an infrared singular ghost dressing function, 
as in Coulomb gauge, while the Type II scenario corresponds 
to a finite ghost dressing function, as appears to be the case in 



Landau gauge. This would most likely be the case if |(|)pA(^)P 
is narrowly peaked around k ~ p, where (^pA (k) is the Fourier 
transform of an F-P eigenmode (^pA (x) with a low-lying eigen- 
value XpA- 

Since the FP spectra at the Gribov horizon have been de- 
rived here from ordinary 2nd order perturbation theory (plus 
an ansatz for the gluon propagator), there is obviously a ques- 
tion of whether perturbation theory can be trusted in this con- 
text. In D = 3 spacetime dimensions the coupling has 
units of mass, so the expansion parameter at p — > will be 
g^/m, while the expansion parameter at large p will be g'^/\p\- 
The perturbative calculation of the FP eigenvalue spectrum at 
p ^ should therefore be trustworthy for large m/g^. Un- 
fortunately, we have seen that the critical mass parameter nic 
corresponding to the Gribov horizon is actually rather small, 
in units of g^, with, e.g., md g^ = 0.223 in Coulomb gauge, 
and mc/g^ = 0.0864 in Landau gauge in three spacetime di- 
mensions and a = 1 . Of course, the perturbative expansion 
may also involve some numerical factors, and without calcu- 
lating to higher orders, or estimating the radius of convergence 
in some way, it is difficult to judge the accuracy of the second- 
order term in the series. But there is no particular reason for 
confidence in the second-order results at m — nic at the quan- 
titative level. It was argued however in section on rather 
general grounds, that it is natural to expect either Type I or 
Type II behavior of the Faddeev-Popov spectrum at the Gri- 
bov horizon. The perturbative calculation, at this stage, sim- 
ply provides a concrete illustration in support of this rather 
general qualitative argument. 
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Appendix 

In order to derive the inequality i2.9i stated in section [III 
we begin with the following 

Theorem. Let Q be any Hermitian operator with a dis- 
crete set of eigenstates {|n)}, whose corresponding eigen- 
value spectrum {qn} is bounded from below, and ordered such 
that qn < qn+i- Here the index n runs from 1 up to the di- 
mension of the Hilbert space Nh (which need not be finite). 
Let {\<^n)} be any other complete set of orthonormal states 
spanning the same Hilbert space as the {|n)}. Then, for any 
N<Nh, 



£((^„IGI^„)>E?« 



(A.l) 



This is a fairly trivial generalization of the inequality under- 
lying the Rayleigh-Ritz variational method (the case = 1), 
and the proof goes as follows: Define 

T = f (^„iei^„> 



EEE(^«l^)(^lGk«)('«i^«) 

«=1 *: m 

Y,qmPN{m) (A.2) 



where 



Observe that 



n=l 



0<PivH<P/v„H = l (A.4) 



and 



Y,PN{m)^N (A.5) 

m= 1 

Since PN{m) < 1, and with regard to the constraint ( lA.SI l, the 
smallest possible value of T is obviously obtained for 



PnH 



1 m<N 
m>N 



(A.6) 



Substituting this optimal choice into the last line of (IA.2b . we 
find that 



T > T 



N 

mm qn 

m=\ 



(A.l) 



and the inequality stated in the theorem is established. 
From this theorem it follows that 

(N) (N) 



where now the {(])„} are the eigenstates of the F-P operator, 
and where we have defined 



{N) 



E^Ee(^'-n-n) 



(A.9) 



n bA 



Now let 



piv(A,) = 



1 



(N) 



NcolorsV „ 



E8(X-X„) (A. 10) 



Assuming non-degeneracy, each X determines n,A uniquely. 



n determines A,!"^, and = A,},"' in the large-volume limit. In 



(0) 



this limit we can therefore we can express (IA.81 ) as 

d'kM'k){if^\^y^\ify}> ldXpNii)pH'>^) (A.ii) 

At small X there is some leading power behavior 



Pn{^)pI) ~ bX" (A.12) 



Then, since (lA.llb must be true for any mode cutoff A^, no 
matter how small, it follows that either r < q, or r ^ q and 
a> b. Either way. 

Finally, given that all the near-zero modes are included in the 
sum ( |A.9l l, we have 



Pn{X)=p(X) as X^O (A.14) 
and the inequality dA.Ob becomes 

lim(Pm^^ (A.15) 

This establishes eq. ( 12. 9b . 
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